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Abstract 



We present the results of a search for pair production of a fourth- 
generation charge -1/3 quark (&') in ^fs = 1.8 TeV pp collisions using 88 pb _1 
of data obtained with the Collider Detector at Fermilab. We assume that 
both quarks decay via the flavor-changing neutral current process b' — > bZ° 
and that the b' mass is greater than mz + m^. We studied the decay mode 
b'b' — ► Z°Z°bb where one Z° decays into e + e~ or [i + n~ and the other decays 
hadronically, giving a signature of two leptons plus jets. An upper limit on 
the Cpp-^wi' x [BR (b' — > bZ )] 2 is established as a function of the 5' mass. 
We exclude at 95% confidence level a b' quark with mass between 100 and 
199 GeV/c 2 for BR(b' -► = 100%. 

13.85.Rm, 13.85.Qk, 14.65.-q 
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The Standard Model (SM) with three generations of quarks and leptons is in excellent 
agreement with all experimental data available today. There is no strong reason to believe 
that an extra fermion generation exists. However, the SM does not explain either the fermion 
family replication or the fermion mass hierarchy. Several models have been proposed to solve 
shortcomings in the SM through the introduction of extra quarks and leptons, while grand 
unified theories, supersymmetry, supergravity and superstrings predict or can accomodate 
extra fermion states An extensive discussion of such models can be found in a recent 
review 0. 

In general, flavor-changing neutral current (FCNC) processes in the Standard Model 
are highly suppressed. However, if a fourth-generation charge -1/3 quark (£/) exists and is 
lighter than both the t' (its partner in an SU(2) doublet) and the top quark (t), the charged- 
current (CC) decays b' — > tW~ and b' — > t'W~ are kinematically forbidden. The leading 
charged- current decay mode will then be b' — > cW~ , which is doubly Cabibbo-suppressed. 
In this situation loop-induced FCNC decays can dominate provided |Kb'|/|^t6'| is l ess 
than roughly 10~ 2 to 10~ 3 , depending on the b' and t' masses jp. If my > mz + mj, the 
dominant FCNC decay mode is b' — > bZ° as long b' — > bH is kinematically suppressed or 
forbidden 0. For m t < my < m t + mw, the decay mode b' — > tW* becomes available but 
is suppressed by three-body phase space, and the b' — > bZ° channel can still dominate over 
the CC decay for b' masses up to about 230 GeV/c 2 

Several experiments have searched explicitly for V quarks decaying via FCNC |7j]. The 
most stringent limit comes from the D0 Collaboration, which searched in the b'b' — > ^gbb 
and b'b' — > ^^bb channels, excluding a b' quark mass up to mz + m^ for a FCNC branching 
fraction larger than 50% M. CDF has excluded a long-lived b' quark with mass up to 
148 GeV/c 2 and a lifetime of r « 3.3 x lO" 11 sec, assuming BR{b' -> bZ°) = 100% [|]. If 
the CC decay b' — > cW~ dominates, the lower mass bound of 128 GeV found in a DO top 
quark search [JTO]] also applies to the b' quark fLl~| . 

In this Letter, we report on a search for a b' quark using 88 ± 4 pb" 1 of pp collisions 
at y^s = 1.8 TeV collected with the CDF detector from 1994 to 1995. Fourth-generation 
V quarks can be pair-produced in pp collisions through gg fusion and qq annihilation with 
the same cross section, for a given mass, as the top quark. We search for pair-produced b' 
quarks decaying via FCNC into bZ°, where one Z° decays into leptons and the other decays 
hadronically. The signature is two high transverse momentum (pt) leptons from the Z° 
decay, two high-p^ jets from the second Z°, and two b jets whose pr scales with the b' mass. 

A detailed description of the CDF detector can be found elsewhere |T2]]. We briefly 



describe the components most relevant for this analysis. Inside a 1.4 T solenoidal magnetic 
field, the silicon vertex detector (SVX), the vertex time projection chamber (VTX), and 
the central tracking chamber (CTC) provide tracking information. The SVX, positioned 
immediately outside the beampipe and inside the VTX, consists of four layers of silicon 



micro-strip detectors and covers \z\ < 25 cm |I3" |. It provides for precise track reconstruction 



in the plane transverse to the beam and is used to identify secondary vertices from the decay 
of b hadrons. The VTX is used to measure the position of the primary interaction vertex 
along the z axis. The CTC is a cylindrical drift chamber that covers the pseudorapidity 
range < 1.1 and consists of 84 layers that are grouped in nine alternating superlayers of 
axial and stereo wires. Outside of the solenoid, electromagnetic and hadronic calorimeters, 
arranged in a projective tower geometry, surround the tracking volume and are used to 
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identify electrons and jets over the range \rj\ < 4.2. The electron energy is measured in 
the central electromagnetic calorimeter (CEM) (\r)\ < 1.1) and the end-plug electromagnetic 
calorimeter (PEM) (1.1 < \r]\ < 2.4). Outside the calorimeters, three systems of drift 
chambers in the region \rj\ < 1.0 provide muon identification. 

We select events satisfying a high-pr lepton trigger, containing a well-identified muon 
or electron in the central region, whose primary vertex is within 60 cm of the nominal 
interaction position. A trigger that requires one jet with Et > 10 GeV, in addition to the 
lepton, is also used for muon events. Inclusive Z° — ► e + e~ and Z° —>■ /x + /x~ samples are 
selected by requiring one primary lepton that satisfies tight lepton identification cuts and 



a second lepton satisfying loose identification cuts ||14|| . Dielectron events are selected by 
requiring at least one tight electron with transverse energy Et > 20 GeV in the CEM and a 
second loose electron with Et > 10 GeV in either the CEM or PEM calorimeters. Dimuon 
events are required to have one tight muon with transverse momentum pt > 20 GeV/c in 
the central region and a second loose muon with > 10 GeV/c. A calorimeter isolation cut 
is imposed on the second lepton. We accept events if the reconstructed ee or /x/x invariant 
mass is between 75 and 105 GeV/c 2 . After this selection there are 6287 (2940) Z° events 
remaining in the electron (muon) data sample. 

In order to optimize our sensitivity to a b' quark signal we make a jet selection that 
depends on the b' mass being considered. Hadronic jets are selected using a clustering 
algorithm |T5| with a cone size of AR = y/Arj 2 + A<p 2 = 0.4. Each event is required to have 
at least three jets within \r)\ < 2.0, two of which with Et > 15 GeV. For b' masses above 
120 GeV/c 2 , the third jet is required to have E T > 15 GeV. For my < 120 GeV/c 2 , the E T 
requirement on the third jet is relaxed to Et > 7 GeV since the b jets for b' masses near 
the rriz + threshold have low momentum. We define the variable E T ets as the summed 
transverse energy of jets with Et > 15 GeV and \r)\ < 2.0 and require this quantity to be 
larger than myc 2 — 60 GeV. Figure [l| shows the J2E T ets distribution for e + e~ and /x + /x _ 
events passing the 3-jet requirement for b' masses above 120 GeV/c 2 . Also shown are the 
distributions expected from SM background and from a V quark with a mass of 150 GeV/c 2 
(see below). 

We further require at least one jet to be tagged as a b quark by the SVX ^-tagging 
algorithm developed for the top quark analysis |16 . The number of events passing each 



major selection criterion for each leptonic channel is shown in Table Q. One /x/x event passes 
all our selection criteria for my < 120 GeV/c 2 . This event has a third jet with Et = 8.3 GeV 
which fails the third-jet Et requirement for larger b' masses. 

The signal acceptance and detection efficiencies are estimated from a combination of 
data and Monte Carlo simulation. We have generated b'V — > bZ°bZ° Monte Carlo samples 
for different b' masses between 100 and 210 GeV/c 2 using the HERWIG program |T7j with 
MRSD0' structure functions ||18|| . One Z° is required to decay into muons or electrons while 
the other is allowed to decay through any available decay channel. The CLEO QQ Monte 



Carlo program |T9[ is used to model the decays of b hadrons. These events are passed 
through a simulation of the CDF detector and subjected to the same selection requirements 
as the data. 

The electron trigger efficiency is determined from data to be (92 ±1)%, while the muon 
trigger efficiency per event (82 ±4)% is obtained from a combination of data and simulation. 
The efficiencies of the lepton identification cuts are determined using a Z° — >• e + e _ (/x + /i~) 
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data sample with an unbiased selection on one of the leptons. The Z° — > e + e~ and Z° — > 
yU+yU" geometric and kinematic acceptance and detection efficiency, including the isolation 
efficiency, is (41 ± 3%) and (30 ±3)% respectively and is nearly independent of the b' mass. 

The event 6-tag efficiency rises with my from 17% for my = 100 GeV/c 2 , to values 
between 50% and 57% for masses above 150 GeV/c 2 . The total acceptance times efficiency, 
not including the BR(Z — ► / + /~), increases from 1.7% (1.6%) to 14% (11%) for the electron 
(muon) channel as my increases from 100 to 210 GeV/c 2 (Table |TJ). This increase is due 
to the fact that a more massive V leads to a more central event with more energetic jets in 
which, in addition, the 6-tag algorithm is more efficient. 

The dominant systematic uncertainties on the acceptance times efficiency arise from 
the jet energy scale and gluon radiation ]14|]. By varying parameters in the Monte Carlo 
simulation we estimate that the systematic uncertainty due to the jet energy scale in the 
electron (muon) channel is 16% (14%) for my = 100 GeV/c 2 and less than 13% for higher 
masses. The presence of gluon radiation increases the jet multiplicity and therefore increases 
the efficiency of the three-jet requirement. This effect is more pronounced at low h' mass 
because the b quarks from low-mass b' decay are produced near threshold and therefore are 
detected with low efficiency. We estimate the systematic uncertainty due to this effect to 
be 19% (18%) in the electron (muon) channel for my = 100 GeV/c 2 , and less than 9% for a 
heavier b'. Other important systematic uncertainties arise from the 6-tag efficiency (10%), 
parton distribution function (5%), total integrated luminosity (4.1%), lepton identification 
efficiency (4% for electrons, 5% for muons), isolation efficiency (4%) and trigger efficiency 
(1% for electrons, 5% for muons). The total uncertainty on the acceptance times efficiency 
is shown function of b' mass in Table |IT| 

The only non-negligible background is from Z° events with associated QCD hadronic 
jets. This background is estimated using a combination of the VECBOS [fZIJ and HERWIG 
Monte Carlo programs. VECBOS is used to make the calculation of the leading-order 
matrix elements for Z° + three partons events, with the MRSD0' structure functions and 



the (pr) of the generated partons for the QCD renormalization and factorization scales [21 



A partial higher-order correction to the tree-level diagrams is obtained by including gluon 
radiation and hadronic fragmentation using HERWIG. These Z° events are then passed 
through a simulation of the CDF detector. We estimate the 6-tag rate in Z° plus jet events 
directly from data using a technique developed for the top analysis f2"2]| . We apply the 



6-tag rates measured in an inclusive jet sample to the Z°+ jets events that pass all the 
other selection criteria. This method overestimates the background because the inclusive 
jet sample contains heavy-quark contributions that are not present in Z° + jets events. We 
expect approximately two background events for my < 120 GeV/c 2 and less than one event 
for my > 120 GeV/c 2 , in agreement with the number of events observed in the data. 

Under the assumption that the observed event is from signal, that is, without sub- 
tracting background, we obtain a conservative 95% confidence level upper limit on the 
a pp~*b'b' x [BR{V — ► bZ )] 2 . The limit is presented as a function of the b' mass in Table ||. 
We have used a Bayesian method to calculate the limit and treat the number of expected 
signal events as a Poisson distribution convoluted with a Gaussian systematic uncertainty. 
Using the theoretical next-to-leading-order b' pair production cross section [2^] and assuming 
that BR{b' — > bZ°) is 100%, we exclude at 95% confidence level b' masses from 100 GeV/c 2 
to 199 GeV/c 2 , as shown in Figure |2|. This search is also sensitive to other b' decay channels 
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such as b' — > bH or b' — > clU~ as long BR(b' — > 6Z) is not negligible, since the hadronic 
decays of the H or W are kinematically similar to those of the Z. The acceptance for 
b'b' — * bbZH is 1.7 to 0.5 times the acceptance for b'b' — > bbZZ, depending on the Higgs 
and b' masses and not including the BR(Z — > Z + Z~). However, if we conservatively assume 
no sensitivity to these decay modes, we exclude a b' mass from 104 GeV to 152 GeV for 
BR(b' -> 6Z) > 50%. 
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FIGURES 




FIG. 1. Ej? ts distribution for events with at least 3 jets with Et > 15 GeV and \t]\ < 2, 
before the 6-tagging requirement. The expected SM background is shown shaded. The expected 
signal event distribution for a b' quark mass of 150 GeV/c 2 is shown as a solid line. The vertical 
dashed line represents the ^2 Ej^ ts cut for this specific b' mass. Events to the right of this line are 
accepted. 



10 



-Q 

NT 
-Q 

T 

V 

DC 
QQ 

D 



10 - 



t — I — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



_L 



_L 



■- Theory a (pp^b'b'X) 



95% CL limit 



199 GeV/f 



_L 



100 120 140 160 180 200 220 

b' mass (GeV/c 2 ) 



FIG. 2. The 95% confidence level upper limit on pp — > b'b'X production cross sec- 
tion times the b' — > bZ° branching ratio squared (solid). The dashed curve shows the pre- 
[BR(b' -> bZ )] with the NLO production cross section from Ref. 23 and 



dieted a, 



pp— ►&'&' 



./ x 



BR(b' bZ°) = 1. 
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TABLES 



TABLE I. Events observed in data after each main selection requirement in both the electron 



and the muon channels. 



my 




Z° e+e- 






Z° - fi+fi- 




(GeV/c 2 ) 


3 jets 


E E J T ets 


6-tag 


3 jets 


J2E J T ets 


6-tag 


100 


34 


31 





32 


29 


1 


120 


34 


20 





32 


21 


1 


140 


9 


8 





8 


5 





160 


9 


4 





8 


4 





180 


9 


1 





8 


3 





200 


9 


1 





8 


2 






TABLE II. Total acceptance (A) times efficiency (e) and relative systematic uncertainties 
(Stotal) in the electron and muon channels, 95% C.L. upper limit on the pair-production cross 
section times the branching ratio of b' — ► bZ° squared, and theoretical pair-production cross sec- 



tion (23]. 





Z°- 


-> e + e 


Z°- 


-> /i+fi 






my 


(A-e) 


Stotal 


(A-e) 


^total 


a ■ BRl 5%CL 


^theory 


(GeV/c 2 ) 


(%) 


(%) 


(%) 


(%) 


(Pb) 


(Pb) 


100 


1.7 


29 


1.6 


27 


37 


102 


110 


4.6 


21 


4.2 


21 


11 


61.6 


120 


7.6 


20 


6.4 


19 


6.5 


38.9 


130 


8.2 


19 


6.8 


19 


3.8 


25.4 


140 


9.9 


19 


8.3 


19 


3.1 


16.9 


150 


11 


19 


9.2 


19 


2.8 


11.7 


160 


12 


19 


9.8 


19 


2.6 


8.16 


170 


12 


19 


10 


19 


2.5 


5.83 


180 


13 


19 


10 


19 


2.4 


4.21 


190 


13 


19 


11 


19 


2.4 


3.06 


200 


13 


19 


11 


19 


2.4 


2.26 


210 


14 


19 


11 


19 


2.3 


1.68 
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